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ABSTRACT 

Aims. Studying molecular gas properties in merging galaxies gives important clues to the onset and evolution of interaction-triggered 
starbursts. The line intensity ratio can be used as a tracer of how dynamics and star formation processes impact the gas properties. 
The Medusa (NGC 4194) is particularly interesting to study since its jE£ ra tio rivals that of ultraluminous galaxies (ULIRGs), despite 
i-C ■ the comparatively modest luminosity, indicating an exceptionally high star formation efficiency (SFE) in the Medusa merger. 

Oh' Methods. High resolution OVRO (Owens Valley Radio Observatory) observations of the 13 CO 1-0 have been obtained and compared 

with matched resolution OVRO 12 CO 1-0 data to investigate the molecular gas cloud properties in the Medusa merger. 
Results. Interferometric observations of 12 CO and 13 CO 1-0 in the Medusa (NGC 4194) merger show the 1-0 intensity ratio (7?) 
increases from normal, quiescent values (7-10) in the outer parts (r > 2 kpc) of the galaxy to high (16 to >40) values in the central (r < 
1 kpc) starburst region. In the central two kpc there is an east-west gradient in K where the line ratio changes by more than a factor of 
three over 5" (945 pc). The integrated 13 CO emission peaks in the north-western starburst region while the central 12 CO emission is 
strongly associated with the prominent crossing dust-lane. 

Conclusions. We discuss the central east-west gradient in K in the context of gas properties in the starburst and the central dust lane. 
We suggest that the central gradient in R is mainly caused by diffuse gas in the dust lane. In this scenario, the actual molecular mass 
distribution is better traced by the 13 CO 1-0 emission than the 12 CO. The possibilities of temperature and abundance gradients are 
0^ ■ also discussed. We compare the central gas properties of the Medusa to those of other minor mergers and suggest that the extreme 

' and transient phase of the Medusa star formation activity has similar traits to those of high-redshift galaxies. 
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1. Introduction Both high kinetic temperatures and large turbulent line 

widths will decrease the optical depth (r), of the 12 CO and 13 CO 

Studying the molecular cloud properties in starburst, active and 1-0 lines and thus elevate the line ratio <R. Therefore, informa- 

interacting galaxies is important in order to understand the feed- tion on the spatial variation of 9? can be used to identify regions 

■ back mechanisms between star formation, dynamics and the in- of extreme or unusual physical conditions in the molecular gas - 

terstellar medium. Important probes of cloud properties include and global values can help classify galaxies. Aalto et al. (1995) 

molecules such as HCN, HCO + and HNC (e.g. Krips et al. 2008; suggested some general diagnostics of the cloud conditions and 

Gracia-Carpio et al. 2008; Imanishi et al. 2004; Gao & Solomon environment based on global values of H: 
2004; Aalto et al. 2002) that trace the dense {nl 10 4 cm 4 ) star 

forming phase of the molecular gas. _ Small ratioS5 ft ~ 6 are an indication of a normal Galactic 

Alternatively, to study the bulk properties of the molecular disc population of clouds dominated by cool, self-gravitating 

gas one can use the ratio between 12 CO and its isotopomer 13 CO. Giant Molecular Clouds (GMCs). 

Globally, there is a correlation between the 12 CO/ 13 CO 1-0 ratio - Intermediate ratios \Q''R < A5 are associated with the inner 

CR) and the FIR /(60yum)//(100^m) flux ratio (e.g. Young & regions of normal starburst galaxies where the gas is warmer 

Sanders 1986; Aalto et al. 1995). The extreme values of % (i.e. and denser than in the disc. 

•R>20) generally occur in luminous merging galaxies with large - The extreme values H > 20 seem to originate in warm, tur- 

dust temperatures (e.g. Aalto et al. 1995; Glenn & Hunter 2001). bulent, high pressure gas in the centres of luminous mergers 

Within galaxies there is a general trend of increasing 7? towards with highly compact molecular distributions and gas surface 

the central region where the gas is warmer and denser (e.g. Wall densities in excess of 10 4 M pc~ 2 — two orders of magnitude 

et al. 1993; Aalto et al. 1995; Paglione et al. 2001). higher than in typical Milky Way GMCs. Large surface den- 
sities require high pressures in hydrostatic equilibrium and 

low density (n < 10 3 cm 4 ) - I2 CO emitting - gas must be 

Send offprint requests to: S. Aalto supported by large turbulent line widths (P oc n(6V) 2 ). Thus, 
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tco can be reduced to moderate (« 1) values, resulting in 
large H. 

More recent studies show that high values of "R also occur more 
locally in less extreme galaxies. Tosaki et al. (2002) find this in 
the spiral arms of M 5 1 and Meier & Turner (2004) in the centre 
of the spiral galaxy NGC 6946. Huttemeister et al. (2000) find 
large values of "R in the large scale bar of NGC 7479 as a re- 
sult of the dynamical impact of a density wave on gas properties. 

The nearby merger NGC 4194 - the Medusa merger - be- 
longs to a class of lower luminosity (Lfir = 8.5 x 10 10 L Q at 
D=39 Mpc) E+S mergers (Aalto & Huttemeister 2000; Manthey 
et al. 2008), an order of magnitude fainter than well known Ultra 
Luminous Galaxies (ULIRGs (L F ir; 10 12 L )) such as Arp 220. 
The Medusa has an elevated value of % («20) (e.g. Aalto et al. 
1991; Casoli et al. 1992; Glenn & Hunter 2001), but in con- 
trast to the more luminous high-'R mergers, NGC 4194 has a 
relatively extended (2 kpc) starburst region (e.g. Armus et al. 
1990; Prestwich et al. 1994) and Aalto & Huttemeister (2000) 
found that the molecular gas was also distributed on the compar- 
atively large scale of 25" (4.7 kpc), despite its advanced stage 
of merger. Due to the linear extent of the molecular gas and rel- 
ative proximity of the Medusa merger, it is possible to spatially 
resolve the relative distributions of 12 CO and 13 CO and hence 
allow the investigation of the underlying causes of the elevated 
H values observed in this galaxy. The Medusa starburst is par- 
ticularly interesting to study since its ^ ratio of 163 rivals that 
of ULIRGs (for example Arp 220 with a of 210), despite 
the comparatively modest luminosity, indicating an exception- 
ally high star formation efficiency (SFE) in the Medusa merger. 
Gao & Solomon (2004) find that for galaxies with L FTR ; 10 11 L G 
the ratio is typically 33. At higher luminosity they find that 
^ instead increases to 100-300. Gao & Solomon (2004) sug- 
gest that, in contrast to CO, the ratio is independent of lu- 
minosity with |^=900. Interestingly the Medusa deviates sig- 
nificantly from this HCN correlation with a ; 4075 based on 
upper limits from Aalto & Huttemeister (2000)). Further discus- 
sion of the underlying causes of this discrepancy are presented 
in Aalto & Huttemeister (2000). 

We have imaged NGC 4194 in the 13 CO 1-0 line with the 
Owens Valley Radio Observatory (OVRO) millimetre array. Our 
aim was to compare the distribution of the 12 CO and 13 CO emis- 
sion within the galaxy and to see how it relates to star formation, 
gas surface density and dynamics. In sections 2 and 3 we discuss 
the observations and results. In section 4 we discuss the "R line 
ratio variations within the Medusa merger in terms of warm gas 
in the starburst region and possible diffuse gas in the central dust 
lane. Future observational tests for these proposed scenarios are 
suggested. 




Fig.l. The 110 GHz continuum of NGC 4194 overlaid on an HST 
WFPC2 image. The contours are: 0.9 x (-3, 3, 4, 5) mjy beam -1 . Rms in 
the map is 1 mjy beam -1 thus the first contour is at 2.7cr. The peak flux 
density is 4.6 mJy beam -1 . The cross marks the position of the 1 .4 GHz 
continuum (Condon et al. 1990). Hatched ellipse indicate beam size and 
orientation. 



1 1 70 km s - 1 . Data were binned to 4 MHz resolution, correspond- 
ing to 1 1 km s -1 , and to construct the map, the resolution was re- 
duced to 44 km s -1 . The sensitivity of this map is 3 mJy beam" 1 
channel -1 , corresponding to 0.014 K channel -1 . 



2.1. Continuum subtraction 

Using the aips task imlin we fitted continuum to the line-free 
channels and subtracted from the cleaned map. The resulting 
continuum map is presented in Fig. 1 and the flux density and 
position in Tab. 1 . The 1 .4 GHz continuum position of Condon 
et al. (1990) agree within 0."4 with the 1 10 GHz position which 
is within the positional errors of the 13 CO map. The 1.4 GHz ra- 
dio continuum emission is dominated by synchrotron emission 
resulting primarily from star-formation, along with a small con- 
tribution from the compact core (Beswick et al. 2005; Condon 
et al. 1990). The 1 10 GHz the emission (see Fig. 1) is coincident 
with the area of highest dust obscuration observed in the optical. 
This emission is tracing the thermal emission from dust in this 
region illuminated by ongoing star-formation. This is consistent 
with the alignment in the 110 and 1.4 GHz continuum positions 
observed. 



2. Observations 

We have obtained maps of 13 CO 1-0 using the Caltech six- 
element OVRO millimeter array. Two tracks were taken in 
the low resolution configuration in April 2000. The naturally 
weighted synthesised beam size is 4'.'56 x 3 '.'98 (472 x 378 pc 
for D=39 Mpc) and the beam position angle (BPA) is -40°. 

System temperatures were 500-600 K. The quasar 
1150+497 was used for phase calibration and Neptune and 
Uranus for absolute flux calibration. The digital correlator, cen- 
tred at V L sr = 2560 km s -1 , provided a total velocity coverage of 



2.2. Alignment of the data 

Both the 12 CO and 13 CO observations were made using the same 
instrument and employing identical gain and passband calibra- 
tors, thus minimising potential misalignments due to instrumen- 
tal effects. Furthermore, the 1 10 GHz continuum extracted from 
the 13 CO data set is positionally coincident with the 1.4 GHz 
continuum (see Sect. 2.1) suggesting that there are no system- 
atic errors in the positions of the 13 CO data set. We also note that 
the shift between 12 CO and 13 CO is partially caused by missing 
13 CO emission in some velocity channels (see Sect. 3.3). 
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Fig. 3. Left: 12 CO contours on a WFPC2 image zoomed in on the central region. Centre: 13 CO contours on WFPC2 image. Both 12 CO and I3 CO 
contours are 40 % of peak flux x (1, 1.5, 2, 2.25). For 13 CO the first contour starts at 2.5cr, for 12 CO the first contour starts at 24 a. Right: Line 
ratio map where the greyscale ranges from light C??=10) to dark (7?=60). Both limits and real values in greys (same scale), and the area where 
there are actual measures overlaid with contours (yellow). Contours are ratios K= 15, 25, 35, 45. Limits are at the 3cr l3 CO detection threshold. 
Please note the small pixel-size (0."5) of the 12 CO and 13 CO maps which may give the impression that there is structure beyond the actual angular 
resolution of the 12 CO and 13 CO maps. 



Table 1. Continuum and I3 CO line results. 



Continuum" : 
position (J2000) 



flux density(mjy) 

Line: 
peak flux density 
(mJy beanT 1 ) 

integrated flux density* 

(Jykms- 1 ) 

(Jykms- 1 ) 



a: 12:14:09.66 
S: 54:31:35.0 

4.6+1 



25 ±3 



3.3+0.6 (centre) 
1.6 + 0.5 (north) 



a) A Gaussian was fitted to the continuum image. The error is lcr rms. 
The 1.4 GHz radio continuum position reported by Condon et al. (1990) 
is: a: 12:14:09.66 6: 54:31:35.5 (J2000). 

b) A Gaussian was fitted to the integrated intensity map of 13 CO. The 
peak integrated flux is 1.9 ± 0.3 (Jy kms~') in the centre. Thus, the 
central source is slightly resolved. The error is lcr rms. 



3. Results: Morphology of the 13 CO emission 

3.1. Previous observations: The n CO distribution 

High resolution OVRO aperture synthesis maps of the 12 CO 1-0 
emission in the "Medusa" galaxy merger (NGC 4194) are pre- 
sented in Aalto & Hiittemeister (2000). It was found that the 
molecular emission is surprisingly extended. The 12 CO emis- 




Fig. 2. Left: Overlay of the l2 CO contours (white, apart from clouds 
in the tail marked in black) on a greyscale, overexposed optical R- 
band image (Mazzarella & Boroson 1993). The dashed curve marks 
the edge of the OVRO primary beam (from Aalto & Hiittemeister 
(2000)). Right: Integrated 13 CO 1-0 line emission in greyscale over- 
laid on a 12 CO contour map smoothed to the 13 CO resolution of 4."5. 
The greyscale range is 0.9 - 2.0 ly km s^'beam -1 . The CO contours are 
5.3 x (1,2,3,4,5,6,7,8) ly km s~'beanr 1 (starting at 6cr). Cross indicate 
position of 1.4 GHz radio continuum peak (as in Fig. 1). 



sion is distributed on a total scale of 25" (4.7 kpc) — despite 
the apparent advanced stage of the merger. The complex, strik- 
ing 12 CO morphology occupies the centre and the north-eastern 
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Fig. 4. Channel maps of 12 CO (left) and 13 CO (right). Boxes indicate channels of "missing" 13 CO emission compared to 12 CO. The 12 CO map 
has been smoothed to the same resolution as 13 CO. Intensity levels for 12 CO and 13 CO are same fraction of peak: 35% (-1.5,-1,1.5,2,2.5) - where 
the peak is 360 and 25 mJy beam-1 for the 12 CO and 13 CO respectively. Note that for 12 CO a big fraction of the detected, extended flux is not 
shown due to the cut-off. The comparative channels maps show the position of the peak emission for 12 CO and 13 CO. The velocity in panel 
18=2714kms~', panel 14 is 2538 kms -1 , and panel 10=2362 kms~'. Increment between panels is 44kms~'. Crosses indicate the 1.4 GHz radio 
continuum position. The beam is marked in the upper right corner of the first panel for the 12 CO and 13 CO channel maps. Greyscale for the 12 CO 
map ranges from 50 to 400mJybeairr' and for 12 CO from 7.5 to 26mJy beanT 1 . For the 13 CO map the first contour corresponds to 2.5<r. 



part of the main optical body (see Fig. 2 (left panel)). The ex- 
tended 12 CO flux density traces the two prominent dust lanes: 
one of which crosses the central region at right angles to the op- 
tical major axis, and the second which curves to the north-east 
and joins the base of the northern molecular tidal tail. The total 
molecular mass in the system is estimated to be at most 2 x 10 9 
M , depending on which 12 CO - H2 conversion factor is appli- 
cable. 



3.2. Integrated 13 CO 1-0 intensity 

The naturally weighted, 4. "5 resolution, 13 CO emission is dis- 
tributed on a scale of 8" (1.8 kpc) in an oval structure centred 
on the optically bright centre of the galaxy. There is also an off- 
centre 13 CO concentration 10" north of the centre. In Fig 2 (right 
panel) we show the integrated 13 CO map overlaid on a 12 CO map 
smoothed to the 13 CO resolution. Integrated- and peak flux den- 
sities are listed in Tab. 1 . 

The 13 CO emission in the inner 8" is shifted towards the 
north-west as compared to the 12 CO distribution. The shift is 
2'.'3 (430 pc). Also the northern 13 CO emission peak is shifted 
off the corresponding 12 CO peak to the south-west by 2'.'5. In 
Fig. 3 (centre panel) the 13 CO 1-0 map is overlaid on an HST 
WFPC2 image. This shows that the 13 CO emission is not tracing 
the central dust lane (which the 12 CO does (left panel in Fig. 3)) 
but instead is shifted towards the north-west part of the starburst 
and towards the western dust features. The moment maps were 
produced through smoothing in velocity and space (3 channels, 
3 pixels boxcar smoothing) and then applying a 2cr flux cut-off 
(allowing both positive and negative flux levels). 



3.3. The channel map 

In Fig. 4 the channel maps of 12 CO and 13 CO 1-0 are pre- 
sented, smoothed to a resolution of 5". Compared to 12 CO, 
13 CO-emission is "missing" primarily in three velocity bins in 
the centre: 2495, 2538 and 2620 km s . In the first case we find 
13 CO emission 10" north of the centre - despite 12 CO in the same 
bin having a peak in the centre. In the second case I3 CO emis- 
sion is found only west of the centre while I2 CO is distributed 
east and west of the radio peak. In the 2620 bin there is a bright 
12 CO peak in the centre which is missing in I3 CO. In general, the 
channel maps show that there is no 13 CO emission east, north- 
east of the radio continuum peak apart from a small region of 
emission 10" north of the central region. 

3.4. Line ratios 

We find that % is 7 - 10 in the region 10" north of the centre - 
while in the centre it is significantly greater: ranging from 16 in 
the western region of the centre to >45 (3cr) in the eastern part 
of the centre where the central dust lane is making a 90 degree 
turn. Over a distance of 5" (945 pc) the line ratio, %, changes 
by more than a factor of 3. A map of the varying fl line ratios is 
displayed in Fig 3 (right panel). 



4. Discussion 

The line ratio change in the Medusa occurs on scales of 0.5 to 1 
kpc where the ratio goes from extreme ratios >40 to more nor- 
mal starburst ratios of 16. This observed gradient in % is likely to 
be caused either by the presence of diffuse, unbound molecular 
gas or a temperature gradient in the ISM. In the following sec- 
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tion we discuss each of these scenarios and compare these with 
observations in other galaxies. 

4.1. Diffuse gas in the central dust lane? 

Optical images of the Medusa merger are dominated by the 
strong absorption from several dust lanes - the most prominent 
one crossing the central starburst region. We suggest that the in- 
tense starburst activity is being fed by gas funnelled along the 
dust lane - supported by the bright 12 CO emission associated 
with the dust lane(s), and the distribution and dynamics seen in 
the Hi absorption at high resolution (Beswick et al. 2005). The 
central dust lane is in front of the main burst region in the east 
and then appears to curve into the burst region itself at the cen- 
tre. A possible scenario is then that the molecular gas is diffuse 
(non self-gravitating) in the dust lane region - creating the gen- 
eral east-west gradient in the line ratio ft. 

Interestingly, a comparison of the Medusa 13 CO image and 
the Ha images of Hattori et al. (2004) show that the position 
of the Ha equivalent width (EW) peak coincides with the re- 
gions of lower values of ft. The large EW is consistent with the 
region being an extra-nuclear star forming region. It coincides 
with the western radio continuum arm and the Giant Molecular 
Association (GMA) feature c in the high-resolution CO map 
(see Aalto & Huttemeister (2000)). Thus, it seems that the 13 CO 
emission correlates spatially with more quiescent, spiral-arm 
like star formation, shifted from the dust lane and the central 
starburst. Note that if the 12 CO 1-0 is mainly tracing diffuse gas 
the 13 CO 1-0 emission is showing the real mass distribution of 
the molecular gas. Fig. 10 of Aalto & Huttemeister (2000) shows 
that the correlation between 1 .4 GHz radio continuum and the 
central dust lane is poor - this could be a further indication of 
both a lack of ongoing star formation in the dust lane as well as 
the molecular gas potentially being diffuse in the dust lane (see 
for instance the discussion of 12 CO and radio continuum in the 
spiral arm of M 83 by Rand et al. (1999)). 

4.1 .1 . Spatial shifts in ft in other galaxies 

Even if the central region of the Medusa is more chaotic than 
that of a density wave spiral galaxy, and the extinction struc- 
ture suggests a three-dimensional central region, it is tempting 
to compare the shift in 12 CO and 13 CO to that found in other 
galaxies. Tosaki et al. (2002) found a spatial shift between 12 CO 
and 13 CO 1-0 emission on similar spatial scales in the southern 
spiral arm of M 5 1 . The shift seems to separate diffuse (unbound) 
lower density gas (as traced by 12 CO) from self-gravitating gas 
emitting the bulk of the 13 CO line emission. Tosaki et al. (2002) 
suggested that there is a 10 7 yr time delay between the accumu- 
lation of gas by the M 5 1 spiral density wave and the formation 
of self-gravitating clouds, resulting in the 13 CO emission being 
found downstream in the spiral arm where it is also spatially 
correlated with Ha emission. Meier & Turner (2004) find for 
NGC 6946 that the value of ft reaches high values of 40 away 
from the central starburst which they attribute to diffuse, low- 
density molecular gas in and behind the molecular arms. 

Within the barred galaxy NGC 7479 spatial shifts between 
13 CO and 12 CO have also been found (Huttemeister et al. 2000). 
These shifts have been attributed to density-wave or bar-like dy- 
namical effects. The similarity of these spatial shifts with those 
found in the Medusa may suggest similar dynamical effects may 
be important. In particular, the transition region between diffuse 
inter arm gas and self-gravitating gas in the arms of M 5 1 sug- 



gesting that the similar shift in the Medusa is caused by the fun- 
nelling and compression of gas towards the starburst region. 

4.2. Temperature gradient in the gas? 

In this scenario, the elevated value of ft in the eastern part of the 
central region is caused not by diffuse molecular gas but by the 
gas there being warmer and denser than in the north-west. For 
a temperature gradient to result in elevated values of ft it is re- 
quired that the average gas densities exceed 3 x 10 3 cirT 3 (since 
the 1-0 transition of 12 CO and 13 CO must be thermalised). A 
change of ft from 16 to >40 can be caused by a temperature in- 
crease from 40 to 150 K in the gas (for a constant column- and 
gas density per cloud). If both the temperature and gas density 
is increasing towards the east then a smaller change in tempera- 
ture could explain the change in ft. If it is filled with warm and 
dense gas, the dust lane would be more likely to harbour intense 
star formation than the rest of the central region. However, as we 
mentioned above, the dust lane shows no sign of elevated star 
formation rates, although it is potentially hidden behind large 
masses of dust. High resolution imaging of higher /-transition 
(e.g. 2-1) emission from 13 CO and 12 CO will allow these diffuse 
gas and temperature effects to be distinguished. If the faint 13 CO 
1-0 emission in the eastern central part is caused by diffuse gas 
the 13 CO 2-1 emission should be even fainter there while the 
opposite is true if it is caused by a high gas temperature. This 
is because the 2-1 level would be more populated than the 1-0 
for 13 CO if the gas is warm - while the opposite is true if the gas 
is diffuse. One example where faint 13 CO 1-0 is accompanied 
by bright 13 CO 2-1 is the merger Arp 299 where it was sug- 
gested that the elevated value of ft is due to temperature effects 
(Aalto et al. 1999). Other studies of multi-transition 13 CO can 
be found in Aalto et al. (e.g. 1995); Glenn & Hunter (e.g. 2001); 
Israel (e.g. 2009). Additionally, high density tracer species such 
as HCN or HCO + will reveal where gas of densities n > 10 4 
crrT 3 is located in relation to the starburst and dust lane. If the 
dust lane is filled with diffuse gas very little dense gas will be 
found there. A high resolution dust SED will reveal where hot 
and cold dust is in relation to the dust lane and central- and off- 
nuclear star formation. 

An abundance gradient: If the physical conditions in the gas are 
kept constant, an elevated ft may be caused by a change in the 
12 C/ 13 C abundance ratio. This possibility has been discussed by 
Casoli et al. (e.g. 1992) as an effect of low metallicity gas being 
transported from the outskirts of the merger to its centre. The 
starburst would then contribute to enriching the gas in 13 C. In 
this scenario the gas in the dust lane has then not yet been en- 
riched by the starburst and there is therefore an east-west and 
north-south age-gradient. However, abundances in NGC 4194 
are found to be close to solar, and no significant metallicity gra- 
dient can be found from a study of the properties of the star form- 
ing regions (Hancock et al. 2006). 

4.3. The large scale gradient in ft 

In the Medusa, the lowest ratios of 7 - 10 are found in the tail- 
like region north of the centre. It is not unusual to find that ft is 
decreasing away from the central region. This has been seen to- 
wards many starburst galaxies including NGC 3256, NGC 1808 
and large spiral galaxies such as M 51. The gas 10" north of the 
centre may therefore consist of rather ordinary cool Galactic- 
type GMCs. This is interesting since the gas there does not ap- 
pear to be in a normal disc distribution, but rather in a shell or 
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a tidal tail. There is no apparent star formation going on in this 
region and the line widths are narrow. If the gas in the central 
region is accumulating from the outer regions of the Medusa via 
the central dust lane, then our results suggest that the molecu- 
lar cloud properties are normal in the outskirts, before becoming 
diffuse due to the funnelling of gas towards the centre and then 
accumulates as self-gravitating, star-forming clouds in the centre 
of the Medusa. 



4.4. The nature of the Medusa starburst 

4.4.1. Comparing with local minor mergers 

One of the few E+S minor mergers to have its molecular medium 
studied is the Medusa "look-alike" NGC 4441. Based on single 
dish observations of 12 CO and 13 CO 2-1 and 1-0 Jiitte et al. 
(2010) suggest that the molecular gas is in an unbound, dif- 
fuse state. Interestingly NGC 4441 appears to be an evolved ver- 
sion of the Medusa with similar morphology and luminosity and 
high resolution 12 COl-0 observations of NGC 4441 show that 
it traces a centrally crossing dust lane. The star formation rate 
of NGC 4441 is much lower than that of the Medusa - likely 
due to the diffuse state of its molecular medium. In this paper 
we suggest that a fraction of the Medusa molecular ISM may be 
the diffuse - it is possible that the entire molecular gas content 
of NGC 4441 is in this state. This raises interesting questions on 
the nature of the Medusa starburst and what it will evolve into. 
Perhaps it will be left with a remnant, feeble dust lane, a small 
reservoir of diffuse gas and only low level star formation, like 
NGC 4441. 



4.4.2. Is the star formation out of equilibrium? 

Despite being a minor merger of comparatively moderate lumi- 
nosity, the Medusa has an L F tr/Lco ratio (163) similar to those 
of typical ULIRGs such as Arp 220 (210) - suggesting that its 
SFE rivals that of the compact ULIRGs. However, the fraction 
of dense (« > 10 4 cirT 3 ) gas is significantly lower in the Medusa 
despite its similar SFE to ULIRGs (Aalto & Hiittemeister 2000). 
Aalto & Hiittemeister (2000) estimate the gas consumption time 
for the Medusa starburst to 40 Myrs and, if a significant frac- 
tion of the molecular medium is diffuse, this time could be even 
shorter. 

The picture becomes even more interesting when one con- 
siders that most of the ongoing star formation in the Medusa 
seems not to be traced by FIR or radio emission. The Ha SFR 
is found to be *46 M yr~' (Hancock et al. 2006), while the 
FIR estimated SFR is 6-7 M yr _1 (Aalto & Hiittemeister 2000). 
(Note, that the SFR estimated from Ha would suggest an even 
more extreme SFE than the already high SFE derived from the 
FIR emission.) It is possible that variations in the initial mass 
function (IMF) can account for some of these differences (e.g. 
Wilkins et al. 2008), although the issue is far from resolved. 
Such an anomaly in the form of a flatter upper IMF could lead 
to enhanced Ha emission. The molecular morphology is also 
quite different from that of 1 .4 GHz radio continuum - which is 
also quite different from that found from other luminous galax- 
ies where the two distributions are usually found to coincide. 
The optical star formation is concentrated in knots that could be 
precursors to globular clusters. These knots have no clear spatial 
relation to either the molecular gas distribution or the continuum. 
This indicates a local deviation from the Kennicutt-Schmidt rela- 
tion (KS) in the Medusa that warrants further study. We suggest 



that the starburst of the Medusa is in an extreme transient phase 
of very high efficiency. 

4.4.3. Comparing the Medusa to high-redshifts extreme 
starbursts 

The Medusa merger has a number of features in common with 
some extreme starbursts at high-z (despite its significantly lower 
luminosity). The starburst regions in high-z systems often seem 
extended (> 3 kpc), gas consumption times are similarly short 
and there is a combination of ordered rotation and merger-driven 
random motions and inflows. More importantly, Bothwell et al. 
(2010) find for three z=3 ULIRGs a significant size difference 
between the CO distribution and star formation tracers. They 
find that this size difference results in the SFEs within systems 
to vary by up to a factor of five. As a consequence to their re- 
sults they conclude that SMGs lie significantly above the KS 
relation, indicating that stars may be formed more efficiently in 
these extreme systems than in other high-z starburst galaxies. 
The Medusa merger is (by far) less luminous than the high-z 
systems studied by Bothwell et al. (2010) but a careful study of 
the extreme starburst of the Medusa merger may lead to insights 
that can be applied to more distant galaxies. 

5. Conclusions 

We have obtained a high resolution map of the Medusa merger, 
NGC 4194, in the 13 CO 1-0 line with the OVRO array. The main 
conclusions we draw from this map are: 

1. The 13 CO 1-0 emission is mainly distributed in the north- 
western part of the inner 8" region. There is also an off- 
centre 13 CO peak in the tail-like distribution to the north, 
where the CO/ 13 CO line ratios are similar to those of galac- 
tic type GMCs 

2. The 12 CO/ 13 CO 1-0 line ratio, % varies on a scale of 950 pc 
from 16 (typical values for starburst galaxies) in the north- 
west to >40 in the eastern part of the central region. 

3. We suggest that this shift between 12 CO and 13 CO emission 
is caused by the presence of diffuse gas in a prominent dust 
lane crossing the eastern part of the centre. A strong density 
wave-like phenomenon may cause shocks (hence the dust 
lane) and the 13 CO is tracing more self-gravitating gas in 
the post-shock region in the starburst. In this scenario, the 
actual molecular mass distribution is better traced by the 
13 CO 1-0 emission than the 12 CO. Alternatively, the line ra- 
tio shift is caused by temperature and density gradients in 
the gas. These two scenarios can be tested through imaging 
of 13 CO 2-1 and high density tracer molecules such as HCN 
andHCO + . 

4. The extreme star formation efficiency and the lack of cor- 
relation between various star formation tracers suggest that 
the inner region of the Medusa merger is in a highly tran- 
sient phase. We predict that it will evolve into a quiescent 
object similar to NGC 4441, but that in its current evolution- 
ary stage it may have some features in common with high 
redshift starburst galaxies. 
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